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ABSTRACT 

We searched for cluster X-ray luminosity and radius evolution using our sample of 201 galaxy clusters detected 
in the 160 deg 2 survey with the ROSAT PSPC (Vikhlinin et al. 1998). With such a large area survey, it is possible, 
for the first time with ROSAT, to test the evolution of luminous clusters, L x > 3 x 10 44 ergs s" 1 in the 0.5-2 keV 
band. We detect a factor of 3-4- deficit of such luminous clusters at z > 0.3 compared to the present. The evolution 
is much weaker or absent at modestly lower luminosities, 1-3 x 10 44 ergs s" 1 . At still lower luminosities, we find 
no evolution from the analysis of the logA^-logS relation. The results in the two upper L x bins are in agreement 
with the Einstein EMSS evolution result (Gioia et al. 1990a, Henry et al. 1992) while being obtained using a 
completely independent cluster sample. The low-L t results are in agreement with other ROSAT surveys (e.g. 
Rosati et al. 1998, Jones et al. 1998). 

We also compare the distribution of core radii of nearby and distant (z > 0.4) luminous (with equivalent temper- 
atures 4-7 keV) clusters, and detect no evolution. The ratio of average core radius for z ~ 0.5 and z < 0. 1 clusters 
is 0.9 ± 0. 1 , and the core radius distributions are remarkably similar. A decrease of cluster sizes incompatible with 
our data is predicted by self-similar evolution models for high-Si universe. 

Subject headings: galaxies: clusters: general — surveys — X-rays: galaxies 



1. INTRODUCTION 

The cluster evolution rate is a strong test of cosmological pa- 
rameters (e.g., White & Rees 1978, Kaiser 1986, Eke, Cole & 
Frenk 1996). It is best to study evolution using X-ray selected 
samples of distant clusters which are much less affected by pro- 
jection than the optically selected samples (van Haarlem et al. 
1997). Of all the interesting cluster parameters such as mass, 
velocity dispersion, and temperature, the X-ray luminosity is 
the most accessible to measurements with present-day instru- 
ments, and most of the earlier studies were focused on evolution 
of the cluster X-ray luminosity function. 

A strong evolution of cluster luminosities at z ~ 0.1 was re- 
ported from the EXOSAT survey (Edge et al. 1990), but was 
later disproved by the ROSAT All-Sky Survey (Ebeling et al. 
1997). At higher redshifts, negative evolution of the cluster 
X-ray luminosity function was first reported by Gioia et al. 
(1990a) using the Einstein Extended Medium Sensitivity Sur- 
vey (EMSS; Gioia et al. 1990b, Stocke et al. 1991). Gioia et 
al. and later Henry et al. (1992) compared the cluster luminos- 
ity functions below and above z = 0.3. They found that while 
the number of the low luminosity clusters does not evolve, 
there is a significant deficit of luminous, L x (0.3-3.5 keV) > 
5 x 10 44 ergs s _1 , clusters at high redshift. 

This EMSS result was questioned recently. Nichol et al. 
(1997) reanalyzed the EMSS cluster sample using ROSAT X- 
ray and new optical observations and argued that the evolu- 
tion reported in the original EMSS papers was not signifi- 
cant. Several groups pursued independent searches for distant 
clusters in archival ROSAT PSPC observations. Collins et al. 
(1997) found that the redshift distribution of 35 clusters de- 
tected in their 17 deg 2 survey is consistent with no evolution. 



This contradicted the earlier claim by Castander et al. (1995) 
of a strong evolution in a similar sample; however, the lat- 
ter authors used an X-ray source detection algorithm not op- 
timized for the cluster search. Jones et al. (1998) presented the 
log N — log S relation for 46 clusters from their 16 deg 2 survey 
and found that this relation is consistent with no evolution of the 
L x < 2 x 10 44 ergs s" 1 (0.5-2 keV band) clusters. Rosati et al. 
(1998) derived cluster luminosity functions up to z ~ 0.8 from 
their sample of 70 clusters detected in a 33 deg 2 survey, and 
found no evolution at low luminosities, L x < 3 x 10 44 ergs s . 
However, none of these ROSAT surveys covers an area large 
enough to probe the evolution of the luminous clusters, and 
their no-evolution claims do not contradict the EMSS results. 

Our 160 deg 2 survey (Vikhlinin et al. 1998, hereafter Pa- 
per I) is the first ROSAT survey comparable with the EMSS 
in sky coverage for distant clusters. We are able to test, and 
confirm, the EMSS evolution results even with the incomplete 
redshift data currently at hand. Eventually, when the spectro- 
scopic work is complete, we will be able to characterize the 
luminosity evolution more accurately. In this Letter, we also 
show that the cluster X-ray core radii do not evolve between 
z ~ 0.5 and now. Throughout the paper, we use definitions /-14 
and L44 for flux and luminosity in the 0.5-2 keV energy band in 
units of 10~ 14 ergs s" 1 cm" 2 and 10 44 ergs s" 1 , respectively. We 
also use Hq = 50 km s" 1 Mpc -1 and qo = 0.5. 



2. CLUSTER SAMPLE 

In Paper I, we presented a catalog of 223 extended X-ray 
sources detected in 646 high Galactic latitude ROSAT PSPC 
observations. For each detected source, we measured the X- 
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ray flux and angular core-radius. We optically confirmed 89% 
of detected sources as clusters of galaxies; 8% are false de- 
tections due to point source confusion and 3% still lack opti- 
cal data. In the high X-ray flux range, which is the focus of 
the present work, 80 out of 82 detected sources are optically 
confirmed clusters. We measured and compiled from the lit- 
erature spectroscopic redshifts for 76 clusters. For the rest of 
the optically confirmed clusters, redshifts are estimated with an 
accuracy of Az s» 0.07 by optical photometry of the brightest 
cluster galaxies. All the X-ray and optical data are presented in 
Paper I. Below we use these data to constrain the evolution of 
cluster luminosities and sizes at high redshift, z ~ 0.5. 

Table 1 

All bright clusters lacking spectroscopic redshifts 



Name 




7 b 


Z range" 


v.min 


2137+0026 


27.8 


0.05 


0.00- 


-0.12 


0.5 


1301+1059 


28.1 


0.30 


0.23- 


-0.34 


0.5 


0841+6422 


29.1 


0.36 


0.29- 


-0.40 


0.5 


1722+4105 


29.4 


0.33 


0.26- 


-0.37 


0.5 


1641+4001 


29.4 


0.51 


0.44- 


-0.55 


0.5 + 


1524+0957 


30.4 


0.11 


0.04- 


-0.15 


0.5 


0259+0013 


32.4 


0.17 


0.10- 


-0.21 


0.5 


2348-3117 


32.5 


0.21 


0.14- 


-0.28 


0.5 


0159+0030 


32.7 


0.26 


0.19- 


-0.30 


0.5 


1515+4346 


34.6 


0.26 


0.19- 


-0.30 


0.5 


0050-0929 


36.6 


0.21 


0.14- 


-0.25 


0.5 


2319+1226 


38.2 


0.25 


0.18- 


-0.29 


0.5 


1146+2854 


39.2 


0.17 


0.10- 


-0.21 


0.5 


1124+4155 


40.1 


0.18 


0.11- 


-0.22 


0.5 


0532^1614 


41.1 


0.10 


0.03- 


-0.17 


0.5 


1013+4933 


45.6 


0.17 


0.10- 


-0.21 


0.4 


1142+2144 


45.9 


0.18 


0.11- 


-0.22 


0.4 


0958+5516 


48.2 


0.20 


0.13- 


-0.24 


0.4 


0237-5224 


64.4 


0.13 


0.06- 


-0.20 


0.4 


1418+2510 


75.6 


0.24 


0.17- 


-0.28 


0.4 


1641+8232 


80.5 


0.26 


0.19- 


-0.30 


0.3 + 


1206-0744 


129.0 


0.12 


0.05- 


-0.16 


0.3 


1630+2434 


179.4 


0.09 


0.02- 


-0.13 


0.3 



a X-ray flux in the 0.5-2 keV band, 10~ 14 ergs s" 1 cm" 2 . 
b Photometric redshift and its 95% confidence range. 

c Minimum redshift required for inclusion in the high-L,, high-z subsample 
(see text); this redshift is defined by the observed flux. Clusters, whose actual 
redshifts can possibly exceed this minimum value, are marked +. 



3. DEFICIT OF LUMINOUS CLUSTERS AT HIGH REDSHIFT 

Cluster evolution is detected in the EMSS sample only for the 
most luminous clusters, L44 > 5 in the 0.3-3.5 keV energy band, 
or L 44 ^ 3 in the 0.5-2 keV band (Gioia et al. 1990a, Henry et 
al. 1992). The lower-luminosity clusters in the EMSS sample 
show little or no evolution. We will search for evolution in 
our sample above this limiting luminosity. Although we cannot 
derive accurate luminosity functions with the presently incom- 
plete spectroscopic data, a sample of luminous high-redshift 
clusters can be selected using the observed X-ray flux and the 
conservative upper bound of their estimated redshift. With such 
a sample, one can test the evolution of the cluster luminosity 
function by comparing the number of detected clusters with the 
prediction of the no evolution model. 

3.1. High-Luminosity and High- Redshift Subsample 

A luminosity L44 = 3 corresponds to observed fluxes, /_i4 = 
77.5, 44.3, and 27.8 at redshifts z = 0.3, 0.4, and 0.5, respec- 
tively. Our high-luminosity, high-redshift subsample is defined 



using these limiting fluxes as follows. The cluster flux must be 
/_i4 > 77.5 and redshift must be z > 0.3, or f-14 > 44.3 and 
z > 0.4, or /_i4 > 27.8 and z > 0.5. For this sample definition, 
the lower limit of the luminosity varies with redshift; the varia- 
tions are, however, limited between L44 = 3 and 6 for z < 0.7. 

There are 48 clusters with fluxes f-14 > 27.8 in our sample. 
Spectroscopic redshifts are available for 25 of them; none of 
these 25 clusters satisfies the selection criteria above. The re- 
maining 23 clusters with photometric redshifts are listed in Ta- 
ble [jj Column (4) in this Table shows the 95% confidence inter- 
val of the photometric redshift. The observed flux corresponds 
to the minimum redshift required by our sample definition (col- 
umn 5). It can be seen that all clusters except 1641+4001 and 
1641+8232 can be confidently excluded from the high-L A , high- 
z subsample. We conclude that at most, only two clusters be- 
long to this subsmaple. 

3.2. Comparison with No-Evolution Model Predictions 

To calculate the expected number of observed clusters, we 
integrated the local luminosity function (Ebeling et al. 1997) in 
the appropriate redshift and luminosity range and accounted for 
the survey solid angle as a function of flux (Paper I). We then 
compare these predictions with the observed number of clus- 
ters in different subsamples. The results are presented in Ta- 
ble H For the no-evolution model, we expect 9.3 clusters in the 
high-L r , high-z subsample, where we observe at most 2. Such a 
deviation is significant at more than 99.5% confidence. For the 
negative evolution observed in the EMSS, we predict that this 
subsample should contain 2-3 clusters, in agreement with the 
observed number. Finally, for qo = 0, the no evolution model 
predicts 8.1 clusters; the observed deficit is still significant in 
this case. 

For a consistency check, we compare the number of z < 0.3 
clusters above the same limiting fluxes, f-14 = 77.5, 44.3, and 
27.8, with the prediction of the no evolution model (Table ||); 
in all flux bins, there is an excellent agreement. 



Table 2 

Comparison with no-evolution model predictions 



Subsample 


Predicted 


Observed 


Prob 


high-L v , high-z 


9.3 


2 


< 0.005 


z < 0.3, f- u > 27.8 


32.5 


39 




z < 0.3, f- 14 > 44.3 


20.4 


21 




z < 0.3, f-u > 77.5 


11.3 


10 




z>0.4, f- u > 13.9 


22.2 


18 


0.22 



We also can compare the observed and predicted number of 
lower-luminosity clusters at high redshift. For that, we use a 
subsample of clusters with fluxes in the range 13.9 < f-u < 
44.3 and z > 0.4. This flux range at z = 0.4 corresponds to the 
luminosity range 1 < L44 < 3, which combines the two lowest 
luminosity bins in the EMSS luminosity function at high red- 
shift. To obtain a conservative lower limit of the number of 
observed clusters, we use the lower bound of photometric red- 
shifts. We find 18 clusters, compared with 22.2 predicted by 
the no-evolution model. That is, the evolution in this interval is 
certainly different from that of the high L x clustrs. 

To summarize, we find a large deficit, by a factor of 3-4, of 
L44 > 3 clusters at high redshift, similarly to the EMSS result. 
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FIG. 1. — Fit to the log N — log 5 relation. The observed logW — logS relation (a) was fit by a family of evolution models which combine luminosity and density 
evolutions as (1 +z) PL and (1 +z) PN , respectively. Solid lines in panel b show 68% and 95% confidence intervals of pjv and pi- The short- and long-dashed lines in 
(b) correspond to no-evolution of the comoving volume emissivity of clusters at z = 0.5, obtained by integration of the luminosity function from to oo (short dash) 
n46 ergs s (long dash). Shaded regions in the right panel correspond to parameters which predict either too many (top, > 6.3) or too few (bottom, 



or from 10 42 to 10 



< 0.36) clusters in the high-L r , high-z sample. The no-evolution model is shown as a dashed line in (a). All allowed models predict virtually the same logA' — logS. 
To illustrate the accuracy of the fit, we show two rejected models, ppi = 0, pt = —1 and pn = 2, pi = (lower and upper dotted lines in the left panel, respectively). 
Error bars show statistical uncertainties of the measured logiV— logS at several representative fluxes. Note that the no-evolution model derived by Rosati et al. (1998) 
differes from that shown here by 10-20% since they integrated only to z = 11 and used a different local luminosity function. 



The evolution rate is smaller at lower luminosities, 1 < L44 < 3, 
again similar to the EMSS and other ROSAT surveys (Rosati 
et al. 1998). At still lower L x , our redshift database is very in- 
complete, but additional constraints on cluster evolution can be 
derived from the log N— log S relation, as described below. 

4. LOG N - LOG 5 RELATION FOR CLUSTERS 

The number of clusters as a function of flux (the logA^-logS 
relation) can provide some constraints on evolution without red- 
shift information. To obtain the evolution constraints from the 
log N — log S, we parameterize the evolution of the luminosity 
function as a combination of pure luminosity and pure number 
density evolution, both as powers of (1+z). With this type of 
evolution, the luminosity function at redshift z can be expressed 
through the local luminosity function Fq(L) as: 

F z (L) = (l+z) p »Fo(L/(l+ Z y L ). (1) 

In this equation, and pi parameterize the rate of the density 
and luminosity evolution, respectively. These parameters can 
be constrained by fitting the observed log N— log S relation. 

To predict the model log.V-log,S relation, we integrated the 
local luminosity function (Ebeling et al. 1997), scaled accord- 
ing to eq. ([j]) in the redshift interval < z < 2 and in the lumi- 
nosity interval 0.01 < L44 < 100. The model logA^-logS was 
normalized by the surface density of clusters expected in the no- 
evolution model above 10~ 12 ergs s" 1 cm" 2 , approximately the 
completeness limit of the Ebeling et al. ROSAT All-Sky Survey 
sample. Finally, we multiplied the model log N— log S relation 
by the solid angle of our survey as a function of flux (Paper I). 

We then found the allowed range of parameters and pi 
using the C-statistic (Cash 1979) calculated in the flux range 



(/-14) from 100 to 4, where the solid angle of our survey still 
exceeds 4 deg 2 . The 68% and 95% confidence region is shown 
by solid lines in the right panel of Fig [[} The allowed combina- 
tion of pjv and pi corresponds approximately to a non-evolving 
comoving volume emissivity. The dotted and dashed lines in 
the right panel correspond to no evolution of the volume emis- 
sivity of all clusters (i.e. defined by Pn + Pl = 0), and clusters 
in the range 10 42 -10 45 ergs s" 1 at z = 0.5, respectively. The no- 
evolution model (pjv = 0, pi = 0) is shown by the dashed line 
in the left panel of Fig [y. All the allowed Pn,Pl parameters 
predict virtually the same log./V-logS relations. On the con- 
trary, the rejected models predict log N — log S relations which 
are markedly different from the data (dotted lines in the left 
panel of Figjl]). 

The observed logAf- log S relation favors no evolution of the 
cluster volume emissivity. This indirectly implies that the lumi- 
nosity function at low L A , which dominates the volume emissiv- 
ity, does not evolve. This, however, does not imply that there is 
no evolution at all. For example, the = 4, p^ = -3 model is al- 
lowed. Since the power-law slope of the luminosity function is 
close to -1 .8 (Ebeling et al. 1997), the number of low luminos- 
ity clusters, L44 < 1, does not evolve in this model, while the 
number of L44 > 3 clusters decreases significantly at z > 0.3. 
This behavior is similar to what we find in §|| 



5. EVOLUTION OF CLUSTER SIZES 

In Paper I, we measured angular core radii by fitting surface 
brightness distributions with the /3-model (e.g., Jones & For- 
man 1984). Here we compare the radius distribution for our 
distant, z > 0.4, clusters and for nearby clusters in Jones & For- 
man (1998) sample. To convert angular radii to proper sizes 
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of distant, z > 0.3, clusters, we used both spectroscopic and 
photometric redshifts. The accuracy of the photometric red- 
shift, ±0.07, is sufficient for this purpose. The correspond- 
ing proper size uncertainty is 15% at z = 0.3, which is smaller 
than the statistical uncertainty of the angular radius measure- 
ments, ~ 20%. Jones & Forman fitted both core-radius and (3 
while we fixed (3 = 0.67 for distant clusters. For consistency, 
we converted core radii from Jones & Forman to (3 = 0.67 using 
Eq. (4) of Paper I. Since Jones & Forman found a correlation 
between core-radius and luminosity, we matched the luminos- 
ity ranges by using only clusters with 1 < L44 < 5 (correspond- 
ing to temperatures 4-7 keV) in both samples. The luminosi- 
ties of distant clusters were computed using both spectroscopic 
and photometric redshifts. The accuracy of photometric red- 
shifts is sufficient because radius is a weak function of lumi- 
nosity. The median redshift of the 25 selected distant clusters 
is Zmed = 0.51. The uncertainty of an individual radius measure- 
ment is <; 30%, which includes photon statistics, background 
modeling, and scatter in /3-parameters (Paper I). This is sig- 
nificantly smaller than the intrinsic width of the derived radius 
distributions. 

The core radius distributions for distant and nearby clus- 
ters are remarkably similar (Fig. ||), especially for go = 0. To 
characterize the radius evolution quantitatively, we use the ra- 
tio of median radii in both samples. The median radius is 
240 ± 14 kpc for nearby clusters, and 210 ± 14 kpc (q = 0.5) 
or 230 ± 16 kpc (qo = 0) for distant clusters. The ratio lies in 
the range 0.95- 1.25. Some models predict a much stronger 
evolution of radii. For example, in Kaiser's (1986) self-similar 
models, the cluster radius grows by a factor of ~ 2 from z = 0.4 
to the present, while hydrodynamic simulations with + A = 1 
(Cen & Ostriker 1994) show a factor of 1.5 growth. However, 
neither these models nor our measurements account for cooling 
flows (e.g. Fabian 1994) which can cause underestimation of 
the core-radius. A comparison of core radii is not meaningful 
if the cooling flow fraction changes with redshift. 



6. CONCLUSIONS 

We present a first ROSAT analysis of the evolution of lumi- 
nous, Lc > 3 x 10 44 ergs s" 1 distant clusters. We find a signif- 
icant, factor of 3-4, decrease in the number of such clusters 
at z > 0.3, confirming the detection of evolution in the EMSS 
(Gioia et al. 1990a, Henry et al. 1992). At lower luminosities, 
1-3 x 10 44 ergs s" 1 , the evolution is undetectable, with a de- 
crease in number by a factor of only 1.3 ± 0.2. This is also con- 
sistent with the EMSS and other ROSAT surveys (e.g. Rosati et 
al. 1998). The absence of evolution of low luminosity clusters 
is also supported by the analysis of the logiV-logS distribution 



from which we find that the cluster volume emissivity, domi- 
nated by low-luminosity objects, does not evolve. The observed 
evolution can be reproduced by a model in which the character- 
istic luminosity decreases with redshift, but the comoving num- 
ber density of clusters increases. Such models arise naturally in 
the hierarchical cluster formation scenario (e.g. Kaiser 1986). 

We compare the distribution of core-radii of distant, z > 0.4 
and nearby clusters. We find that the distribution of core radii 
at z > 0.4 is very similar to that in nearby clusters; the average 
radius has changed at z > 0.4 by a factor of only 0.9 ± 0. 1 . A 
stronger change is expected for hierarchical cluster formation 
in a flat universe (Kaiser 1986, Cen & Ostriker 1994). We also 
note that the assumption of no evolution of cluster sizes has 
been essentially used in flux measurements and area calcula- 
tions in several X-ray surveys (e.g. EMSS, Nichol et al. 1997), 
but is only verified here for the first time. 



-q = 0.5 

q =o 



Core-radius (Mpc) 

FIG. 2. — Core-radius distribution for distant, z > 0.4, clusters, derived from 
our survey (solid and dashed histogram for = 0.5 and <jo = 0, respectively). 
The shaded histogram shows the core-radius distribution for nearby luminous 
clusters from Jones & Forman (1998). The angular resolution limit of our sur- 
vey (15") corresponds to 120 kpc at the median redshift of distant clusters, well 
below the peak of the distribution. 



We thank M. Markevitch for useful comments on the 
manuscript. Financial support was provided by the Smithsonian 
Institution and NAS8-39073 contract. HQ acknowledges sup- 
port from FONDECYT grant 8970009 and the award of Presi- 
dential Chair in Science. 



REFERENCES 



Cash, W. 1979, ApJ, 228, 939 

Castander, F. J. et al. 1995, Nature, 377, 39 

Cen, R. & Ostriker, J. P. 1994, ApJ, 429, L4 

Collins, C. A., Burke, D. J., Romer, A. K., Sharpies, R. M., & Nichol, R. C. 

1997, ApJ, 479, LI 17 
Ebeling, H., Edge, A. C, Fabian, A. C, Allen, S. W., Craford, C. S., & 

Bohringer, H. 1997, ApJ, 479, L101 
Edge, A. C, Stewart, G. C, Fabian, A. C, & Arnaud, K. A. 1990, MNRAS, 

245, 559 

Eke, V R., Cole, S., & Frenk, C. S. 1996, MNRAS, 282, 263 
Fabian, A. C. 1994, ARA&A, 32, 277 

Gioia, I. M., Henry, J. R, Maccacaro, T., Moms, S. L., Stocke, J. T., & Wolter, 
A. 1990a, ApJ, 356, L35 



Gioia, I. M., Maccacaro, T, Schild, R. E., Wolter, A., Stocke, J. T., Morris, S. 

L., & Henry, J. P. 1990b, ApJS, 72, 567 
Henry, J. P., Gioia, I. M., Maccacaro, T., Morris, S. L., Stocke, J. T, & Wolter, 

A. 1992, ApJ, 386, 408 
Jones, C. J. & Forman, W. R. 1984, ApJ, 276, 38 
Jones, C. J. & Forman, W. R. 1998, ApJ, submitted 

Jones, L. R., Scharf, C. A., Ebeling, H., Perlman, E., Wegner, G., Malkan, M., 

& Homer, D. 1998, ApJ, 495, 100 
Kaiser, N. 1986, MNRAS, 222, 323 

Nichol, R. C, Holden, B. P., Romer, A. K., Ulmer, M. P., Burke, D. J., & 

Collins C. A. 1997, ApJ, 481, 644 
Rosati, P., Delia Ceca, R., Norman, C, & Giacconi, R. 1998, ApJ, 492, L21 



CLUSTER EVOLUTION 



5 



Stocke, J. T., Moms, S. L., Gioia, I. M., Maccacaro, T., Schild, R., Wolter, A., 
Fleming, T. A., & Henry, J. P. 1991, ApJS, 76, 813 



van Haarlem, M. P., Frenk, C. S., & White, S. D. M. 1997, MNRAS, 287, 817 
Vikhlinin, A., McNamara, B. R., Fiprman, W., Jones, C, Quintana, H. & 



Hornstrup, A. 1998, ApJ, in press ( astro-ph/9803099 ) 
White, S. D. M. & Rees, M. J. 1978, MNRAS, 183, 34. 



